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Abstract 

The classical (i.e. non-quantum) equilibrium statistical mechanics of a Coulomb fluid 
living on a pseudosphere (an infinite surface of constant negative curvature) is considered. 
The Coulomb fluid occupies a large disk communicating with a reservoir (grand-canonical 
ensemble). The total charge Q on the disk fluctuates. In a macroscopic description, the 
charge correlations near the boundary circle can be described as correlations of a surface 
charge density a. In a macroscopic approach, the variance of Q and the correlation function 
of er are computed; they are universal. These macroscopic results are shown to be valid for 
two solvable microscopic models, in the limit when the microscopic thickness of the surface 
charge density goes to zero. 
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1 INTRODUCTION 



This paper is dedicated to Elliott Lieb on the occasion of his 70th birthday. It is a variation on 
a theme to which Elliott has brought major contributions. [1, 2] 

How statistical mechanics is affected by the curvature of space might be of some interest in 
general relativity, and also is an amusing problem per se. A simple case is a two-dimensional 
system living on a pseudosphere, i.e. a surface of constant negative curvature. Unlike the sphere, 
the pseudosphere has an infinite area, and therefore, on a pseudosphere, one can consider the 
thermodynamic limit of some system while keeping a given curvature. A special feature is that, 
for a large domain, the neighborhood of the boundary has an area of the same order of magnitude 
as the whole area of this domain; this feature makes the approach to the thermodynamic limit 
rather different to what happens in a flat space. 

More specifically, the present paper deals with a two-dimensional classical (i.e. non-quantum) 
Coulomb fluid living on a pseudosphere. This is a system of charged particles interacting by 
Coulomb's law, with this law defined on the pseudosphere, i.e. as the solution of the Poisson 
equation written with the pseudosphere metric. The Coulomb fluid is assumed to be in equi- 
librium and confined in a large disk drawn on the pseudosphere. The grand-canonical ensemble 
is used: the fluid can freely exchange particles with a reservoir. Thus, the total charge Q may 
fluctuate. Furthermore, there are charge correlations which, near the circle boundary of the disk, 
can be described as correlations of a surface charge density a. The aims of the present paper are 
to compute the variance of Q and the two-point correlation function of a. It will be shown that 
these quantities are universal (i.e. independent of the microscopic nature of the fluid). These 
universal results will be checked on two exactly solvable models: the two-component plasma, 
made of two species of particles of opposite signs, and the one-component plasma, made of one 
species of particles in a neutralizing background. 

These problems have already been studied and solved in a flat space. For a finite two- 
dimensional Coulomb fluid in a plane, the total charge Q essentially does not fluctuate [3]. 
Furthermore, in the case of a large disk of radius R centered at the origin, the surface charge 
correlation function [4] is given by the universal expression 

<?<<*m°» t - -s^iwf (L1) 

where (3 is the inverse temperature, a(ip) the surface charge density on the boundary circle at 
the point of polar angle ip, and < . . . > T is a truncated statistical average. These results have 
been checked on exactly solvable models. [3, 5] 

Here, the same problems are considered, now on a pseudosphere. In Section 2, some basic 
properties of the pseudosphere and of Coulomb's law on it are recalled. In Section 3, macroscopic 
electrostatics on a pseudosphere is used for determining the variance of the total charge Q and 
the correlation function of the surface charge density a. In Section 4, the results are checked 
on a solvable model, the two-component plasma at a special temperature. In Section 5, the 
results are checked again on another solvable model, the one-component plasma at a special 
temperature. 

2 PSEUDOSPHERE AND COULOMB'S LAW 

Let us recall a few properties of the surface of constant negative curvature called a pseudosphere. 
Such a surface is a two-dimensional manifold, the entirety of which cannot be embedded in three- 
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dimensional Euclidean space. Its properties are denned by its metric. Several sets of coordinates 
are commonly used. 

The one which renders explicit the resemblance with the sphere is (r, ip) with r £ [0, oo[ and 
<p G ] — 7r, 7r], the metric being 

ds 2 = a 2 (dr 2 + sinh 2 rd99 2 ) (2.1) 

where — 1 /a 2 is the Gaussian curvature (instead of 1/R 2 for a sphere of radius R) . The geodesic 
distance s between two points at (r, <p) and (r', is given by 

cosh(s/a) = coshrcoshr' — sinh r sinh r' cos (y? — <//) (2.2) 

In particular, the geodesic distance of the point (r, p) to the origin is ar. The Laplace-Beltrami 
operator is 

_ 1 / 1 a 9 _L ^\ 

a 2 \ sinh r <9r <9r sinh 2 r dip 2 J 

The set of points at a geodesic distance from the origin less than or equal to R = aro will be 
called a disk of radius R. Its boundary will be called a circle of radius R. Its circumference is 

C = 27rasinhr (2.4) 

and its area is 

^ = 47ra 2 sinh 2 (r /2) (2.5) 

It is remarkable that, for a large radius, both the circumference and the area are proportional 
to expro: the neighborhood of the boundary circle has an area of the same order of magnitude 
as the whole area! 

Another often used set of coordinates is (r, tp) with r/(2a) = tanh(r/2). Then, the metric is 

dg2 = (2 6) 



When these coordinates are used, the whole (infinite) pseudosphere maps on a disk of radius 2a, 
the Poincare disk. 

Finally, here it will be convenient to use also the coordinates (e, p) with 

tanh(r/2) = e _e (2.7) 

Then the metric is 

d S 2 = ^V( de2 + V) (2-8) 
sinh e 

and the Laplace-Beltrami operator has the simple form 

A sinh 2 e / d 2 d 2 \ 

a = ^Gp + vJ (2 - 9) 

The Coulomb potential v(s) at a geodesic distance s from a unit point charge obeys the 
Poisson equation 

Ad(s) = -2tt5^ 2) {s) (2.10) 

where 5^ is the Dirac distribution on the pseudosphere. The solution of (2.10) which vanishes 
at infinity is 

v(s) = -lntanh^- (2.11) 
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3 MACROSCOPIC ELECTROSTATICS, CHARGE FLUCTU- 
ATIONS, SURFACE CHARGE CORRELATIONS 

3.1 Two Problems in Macroscopic Electrostatics 

Here are two problems, the solution of which will be needed in the following. On the pseudo- 
sphere, an ideal conductor fills the disk of radius R = aro centered at the origin. 

Capacitance. The first problem, a very simple one, is: What is the capacitance of this 
disk? If the disk carries a charge Q, this charge uniformly spreads on its circumference and, 
by Newton's theorem (which can be easily shown to be valid on a pseudosphere), generates 
on the whole disk the constant electric potential Qv(R) = — Q lntanh(ro/2). Therefore, the 
capacitance is 

C = "lntanhf ^ 

In the large-disk limit ro — > oo, C ~ exp(ro)/2. 

A Point Charge in the Presence of the Disk. The second problem is: A unit point 
charge is located, outside the disk, at point (r',^' = 0). The disk is grounded (i.e. kept at zero 
potential). What is the electric potential <p(r,(p;r') at some point (r,<p), outside the disk? The 
method of images, which can be used for a flat disk, does not seem to work on a pseudosphere, 
and a Fourier expansion will be used. 

The potential due to the unit point charge alone is (2.11). Expressing this potential in terms 
of cosh(s/a), using (2.2), and expanding as a Fourier series in ip (in the case r < r' which suffices 
for our purpose) gives 

J^L 9 csrnVi e'f p _< ^ 

v(s) = e' + Y, y cosly (t <t<t') (3.2) 

where we have gone from the variables r and r' to the variables e and e' defined by (2.7) and its 
analog for e' . Using (2.9), one easily checks that the Laplacian of each term of (3.2) vanishes. 

The full potential in the presence of the disk is obtained by adding to (3.2) terms of 
zero Laplacian symmetrical in e and e': a term of the form Aqc'c, and terms of the form 
At sinh e'£ sinhe£ cos lip. These terms do vanish when r or r' goes to infinity. The coeffi- 
cients Ai are determined by the condition that the potential vanishes on the disk, i.e. when 
t = tq. The result is 

0(r, <p; t') = e' - — + V 2 sinh e'l ( e~ ei - 4^- sinh el) (3.3) 
eo ^ V sinheo^ / I 

where eo is related to ro by the analog of (2.7). 

3.2 Charge Fluctuations 

The disk of radius R = aro centered at the origin is filled with a Coulomb fluid. It can freely 
exchange charges with a reservoir located at infinity (grand canonical ensemble) . If macroscopic 
electrostatics is applicable, what is the variance of the charge Q carried by the disk? 

The same reasoning as in the flat space case [3], using linear response theory, gives for the 
variance 

P(QY = C (3.4) 
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where here the capacitance C is given by (3.1). This result (3.4) just says that the variation of 
the energy Q 2 /2C has the usual thermal average (l/2)(3~ 1 . In the large-disk limit tq — > oo, 



3.3 Surface Charge Correlations 

In three dimensions, macroscopic electrostatics deals with volume charge densities and surface 
charge densities. In the present case of a two-dimensional system (a disk), the analog of the 
surface charge density actually is a charge per unit length on the boundary circle; we shall 
nevertheless still call it a surface charge density a(ip). If macroscopic electrostatics is applicable, 
what is the two-point correlation function of cr(<^)? 

The same reasoning as in the case of a flat space [4], using linear response theory, gives for 
the two-point correlation function 



(2vra)H drdr' J T=T , ^ 



T — T — TQ 



where </>(t, y;r') is the electric potential (3.3). 

The second derivative in (3.6) can be expressed in closed form in terms of the Jacobi theta [6] 
function 6\. For the sake of dealing only with convergent series, from (3.3) one first computes 
the second derivative for r' > To, i.e. e' < eo: 

( d 2 (j)(T^-T')\ , , ./ 1 ^coshe'/ \ 

Y l ' ' r \ ' = -2sinhe sinhe' — +> — - 1 cos lip \ 3.7) 

V drdr' J T=To \2e ^smhe £ ') 1 > 

For taking the limit of (3.7) when r' — ► To, one substracts from and adds to coshe'^/ sinheo^ a 
term e _ ( eo_e ^. The limit of one of the two resulting series is computed after the summation 
has been performed: 

oo ^ 

lim V e'^-^l cos e<p = 5-^- (3.8) 

e'^eo ^ 4 shr i| 

The other series involves (coshe'f/ sinheo^) — e - ^ 0- " 1 ^. It remains absolutely convergent when 
the limit e' = eo is taken in each term. Thus 



f d^^r'^ . 2 / 1 1 " e -W 

, = ~ 2 sinn eo tttt + 2 > =— » t c 

V 3rdr' / r=T , =7T) ^2e 4sin 2 f ^ 1 - e" 2 ^ 

The sum in (3.9) can be expressed in terms of the Jacobi 6\ function, since [6] 



J-y ,9 l = 7T cot ttv + 4vr V q ,, sin 2^vro (3.10) 

Setting v = (p/{2ir) and q = e~ eo in (3.10), and using its derivative with respect to ip in (3.9) 
gives for the correlation function (3.6) the closed form 



/3<a(^)a(0)) T = ^^sinh 2 e 



1 , 1 d *i(£,e-«; 



e + 7rd^^(^,e-^o) 



(3.H) 
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In the flat-space limit a — > oo, To — > 0, for a fixed value of i? = aro, it can be checked that 
(1.1) is recovered. More interestingly, in the opposite limit of a radius R large compared to the 
"curvature radius" a, i.e. when r — > oo and e — > 0, (3.11) takes a simpler form. Indeed, after 
a Jacobi imaginary transformation [6, 7], the 6\ function can be expressed as the series 



\ 1/2 n=oo 



7H / y> 1 



(3.12) 



If < \ip\ < tt, in the small-eo limit (eo ~ 2e _T °), the first two leading terms of the series (3.12) 
are n = and n = 1. When only these terms are kept, to first order in their ratio exp(— 7r|</?|/eo), 
(3.11) becomes 



-2^ eXPl ' 



e T °7r|(/?| 



(0 < \<p\ < tt) 



(3.13) 



4 TWO-COMPONENT PLASMA ON A PSEUDOSPHERE 

The total charge fluctuation and the surface charge correlation have been obtained under the 
assumption that macroscopic electrostatics is valid. In the large-disk limit R^> a, these results 
(3.5) and (3.13) will now be checked on two solvable microscopic models. Such checks are wel- 
come, because a two-dimensional case when macroscopic electrostatics is not valid, unexpectedly 
at first sight, is known: the charge fluctuations in a short-circuited circular condenser [3]. 

Macroscopic electrostatics uses the concept of surface charge density. Actually, in a mi- 
croscopic model, this "surface density" will have some microscopic thickness, and for macro- 
scopic electrostatics to be valid, it is necessary that this thickness be negligible compared to 
the macroscopic lengths. The microscopic model which will be used in the present section is 
the two-component plasma at a special temperature. Its microscopic scale is characterized by 
a fugacity £, with the dimension (length) -2 . In a disk of radius R on a pseudosphere with 
the "radius of curvature" a, there are two dimensionless parameters involving £: (a 2 and £i? 2 . 
Necessary conditions for macroscopic electrostatics to be valid is that both these parameters be 
large compared to 1. Here, for simplicity, the disk is assumed to be large (R 3> a). 



4.1 Review of the General Formalism 

The two-component plasma is a system of two species of particles, of charges ±1. At the special 
inverse temperature (3 = 2, the model is exactly solvable in different geometries, in particular 
on a pseudosphere [9]. For the sake of completeness, the method of solution is briefly revisited, 
in a form simpler than in the original papers, by a generalization of what has been done in the 
case of a one-component plasma [8]. 

In terms of the coordinates (r, ip) (the Poincare disk representation) , the Coulomb interaction 
(2.11) between two unit point charges at r, and rj is 

v(s) = — In 

where Zj is the complex coordinate of particle j (and zj its complex conjugate): zj = r^e 1 ^'. 
The interaction (4.1) happens to be the Coulomb interaction in a flat disk of radius 2a with 
ideal conductor walls at zero potential. Therefore, one can use the techniques which have been 
developed [10, 11] for dealing with ideal conductor walls. 



(zi - Zj)/(2a) 



1 - 



Aa 2 



(4.1) 



6 



When (3 = 2, the Boltzmann factor for N + positive particles with vector coordinates r+ 
and corresponding complex coordinates zf, 1 < i < N+, and iV_ negative particles with vector 
coordinates r~ and corresponding complex coordinates z~, 1 < i < N_, can be written as (with, 
for the time being, 2a taken as the unit of length: 2a = 1) 

Bn + ,n_ = 

N + N_ 

n (z+- z +)(z+-z+) n -*,-) n n (i-z+ z -)(i-z+ z -) 

l<i<j<N + l<fc<!<AT_ m=ln=l 

JV + AT_ (4-2) 

n (i-z+z+)(i-z+z+) n (i-^5r)( 1 -^ z D n n )(*fc-*n 

l<i<j<iV + l<fe<i<JV_ m=ln=l 

(in the cases iV + = and iV_ = 0, the corresponding products in (4.2) should be replaced by 1; 
in particular B(0, 0) = 1). It is convenient to define 

b n + ,n- - -n^-n: — ( 4 - 3 ) 

Tl Y\( l - z mZm){l - Z n Z n ) 
m=l n=l 

B' is the Boltzmann factor in a disk with ideal conductor walls at zero potential, including now 
in its denominator the contribution from the interaction of each particle with its own image. B' 
has the advantage that it can be written as a N x N determinant (N = N+ + N- is the total 
number of particles), by using the Cauchy identity 

]1 fa -Uj){ Vi -vj) 

_ = (_l)"("-i)/2 det (^_] (4 .4) 
\Ui Vj J i i=1N 



l<i<j<N 



N N 

n n(ui-vj) 
i=ij=i 



Indeed, choosing 

Ui = zf, Vi = l/zf, l<i<N + 
u i+ N+ = l/^ r v i+N+ = z~, 1 < i < iV_ (4.5) 

in (4.4) gives, after some simple manipulations and the reestablishment of an arbitrary value of 
2a, a N x N determinant 

where 



b 'n+,N- = det Aij (4.6) 



^ = ^4zJ {il ^^ N + 

= + 2Q _ if 1 < i < N+, N + <j<N 

z i ~ Z j-N+ 
2a 

= — — if N + < i < N, 1 < j < N + 

Z i-N+ ~ Z j 
4a 2 

Aij = 2 if N + < i,j < N (4.7) 

4a Z i-N+ Z j-N+ 

If the (perhaps different) fugacities are ( + and £_ for the positive and negative particles, respec- 
tively, the grand partition function can be written as 

oo oo .. . -N+ N- 

E = E E JfWl / II d2 *rrMr+) II d 2 r-( + (r-)B' N+!N _ (4.8) 

N+=ON-=0 +' ' J m=0 n=0 
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Indeed, one of the factors [1 — (r 2 /4a 2 )]~ 1 in the area element on the pseudosphere dS = 
[1 — (r 2 /4a 2 )]~ 2 d 2 r has been incorporated into the definition (4.3) of B', while the other factor 
[1 — (r 2 /4a 2 )] -1 has been incorporated in the definition of position-dependent fugacities 



C±(r) 



1 



(4.9) 



Although the integrals in the grand partition function(4.8) diverge (as the separation between a 
positive particle and a negative one goes to zero), this grand partition function can be formally 
manipulated for providing finite correlation functions. 

It will now be shown that the grand canonical partition function can be expressed as one 
determinant of an infinite matrix, continuous in coordinate space. First, one considers the 
functional integral 



Z = J Vm /exp J Y, ^(r)(M-V(ry)yv(r')d 2 rdV 

s,s'=± 



(4.10) 



The fields ip and ip are two-component Grassmann variables (anticommuting variables). The 
components of ip are called ip+ and ip-, and similarly for ip. The covariance of the Gaussian 
measure in (4.10) is the inverse of the kernel M _1 , which is chosen such that 



(Var)Vv(r')> =M ss ,(r,r') 



(4.11) 



where (. . .} denotes an average taken with the Gaussian weight of (4.10) and the 2x2 matrix 
M is 

Second, one considers the functional integral 



4a 2 — zz' 



Z = J VipVijj /exp J Y M^)(M- 1 ) ss ir,r')il; s/ (r')d 2 rd 2 r' 

s,s'=± 



(4.13) 



+ J [(+(r)iP + (r)iP + (r) + C-(r)^_(r)V-(r)]d 2 r 
and one expands Z/Zq in powers of C+( r ) an d C-( r ) as 

u Af + =0Af_=0 + ' ' J m=0 n=0 

x <<Mr+)V> + (r+) • • • Vi + (r+ + )^ + (r+ + )^_(r+)^-(r+) • • • ^(r N J^(r N J) (4.14) 

Third, from the Wick theorem for anticommuting variables [14] and the covariance (4.11), it 
results that the average in (4.14) is equal to the determinant of the matrix A%j defined in (4.7), 
i.e to B' N+ N _ as given by (4.6). Therefore (4.14) is identical to (4.8). The grand partition 
function of the Coulomb gas is 

Z 

E = — (4.15) 
Zq 
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Finally, Z = det(M _1 ) and Z = det(M _1 + Q. In these determinants of infinite order, the 
matrix elements of M are labeled both by the discrete charge indices s, s' and the continuous 
indices r, r'. The infinite diagonal matrix ( is defined as 

Therefore (4.15) does give the grand partition function as the determinant of an infinite matrix, 
continuous in coordinate space: 

Z = det[M(M~ 1 + ()} =det(l + MC) (4.17) 

For computing the densities and many-body densities, some definitions are needed. Let us 
define 

G = (1 + MQ- l M/{A^a) (4.18) 

(the factor [Aira) is there just for keeping the same notation as in previous papers). Thus, G is 
the solution of (1 + MQG = M / (4ira) or, more explicitely, G obeys the integral equation 

G(r, r') + J M(r, r")C(r")G(r", r ')dr" = ^M(r, r') (4.19) 

where it should be remembered that G, M, £ are 2x2 matrices. We also define 

(2 \ 1/2 / 12 \ 1/2 

l -h) &r - r 'A l -h) <4 - 20) 

On (4.12), one sees the symmetries JW ss /(r,r') = ss'M s i s (r', r). By formally expanding the 
definition G = (1 + M()~ l M/(4:7ra) in powers of M( one finds that G has the same symmetries, 
which also hold for G: 

G s Ar,r') = s S 'G s/s (r',r) (4.21) 

The density n s (r) of particles of sign s is given from the grand partition function by a 
functional derivation: 

4^ ) ^ VJ 5( s (r) 

where the factor [1 — (r 2 /4a 2 )] 2 insures that n s (r)dS is the average number of particles in the 
area element dS = [1 - (r 2 /4a 2 )]- 2 d 2 r. Since, from (4.17), lnH = Tr(l + MQ, (4.22), (4.18), 
and (4.20) give 

n s (r) =4vrC s aG ss (r,r) (4.23) 

(actually, for point particles, this density is infinite, but it can be made finite by the introduction 
of a small hard core). The two-body density Ursell functions are given by 

/ r 2 \ 2 / ,2 \ 2 A 2 l n - 

Taking into account the symmetry relations (4.21) gives 

U ss ,(r,r') = - SS '(47rC s a)(47rC s 'a)|G ss ,(r,r')| 2 (4.25) 

From now on, we restrict ourselves to the case of equal fugacities C+ = C- = C- I n t ne Poincare 
disk representation, the Coulomb fluid fills a disk of radius tq. Thus, ((r) = when r > tq. 



n s (r) = (l-^) Cs(r)^k (4.22) 
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The radius r is related to the geodesic radius R = aro by r = 2atanh(ro/2). Without loss of 
generality we can choose the polar angle of r' as <p' = 0. 

The integral equation (4.19) can be tranformed into a differential one, by the application of 
the operator ft = a x d x + cr y d y , where a x and a y are Pauli matrices: 

[ft + AnaC(r)}G(r,r') = 5^ t (r-r') (4.26) 

(2) 

where 5^ t is the Dirac distribution in the plane. This differential equation is to be supplemented 
by the condition that G(r, r') be continuous at the disk boundary r$ and by the boundary 

condition, seen on (4.19), that when r = 2a, G |_ = e 1¥, G++ (and a similar boundary relation 

between (5 H and G ). 

In the case of an infinite system, eq.(4.26) could be solved [9], for r' = 0, in terms of 
hypergeometric functions. In the present case of a finite disk, i.e. when r$ < 2a, an exact 
explicit solution of (4.26) for an arbitrary fugacity seems difficult to obtain. Fortunately, here 
we only need the large-fugacity limit, in which case there are important simplifications. 

4.2 Large Fugacity 

For a flat system, the Coulomb interaction (2.11) becomes — ln(s/2a) where 2a is an irrelevant 
length scale which only contributes an additive constant to the potential. In the flat case [12], 
the rescaled fugacity m = 4ir(a (which has the dimension of an inverse length) was introduced, 
and the correlation length was found to be of the order of vrT 1 . In the present case of a system 
on a pseudosphere, it is convenient to keep the same definition of m. 

On a pseudosphere, in the large fugacity limit AitQa 2 = ma S> 1, if we are interested in 
a solution of (4.26) only in a region of size m _1 , the curvature can be neglected and the flat 
system solutions can be used, with appropriately rescaled coordinates. In particular, if both r 
and r' are sufficiently close to ro, the variation of £(r) can be neglected: in (4.20) and (4.26), 
£(r) can be replaced by the constant C( r o)- Here, we assume the disk to be large, and therefore 
C( r o) ~ C eT °/4- Furthermore (4.20) becomes 

G(r,r') = T G(r,r') (4.27) 

In terms of the rescaled variables (e T °/4)r = t and (e r °/4)r' = t', (4.27) and (4.26) do give the 
flat system equation [12] 

[& +m ]G(r,r') = 4 2 i(t-t') (4.28) 

In an infinite system, the (++) and ( — h) elements of the solution of (4.28) would be 
G ++ (r,r') = (m/2ir)K (m\t - t'|) and G_+(r,r') = (m/27r)e^Ki(m|t - t'|) where K and 
K\ are modified Bessel functions and ip is the argument of te lip — t'. In the present case of a 
finite system in a disk, a "reflected wave" must be added. As a Fourier series in <p, G ++ is of 
the form 

oo 

G ++ (r,r') = ^ [h{mt')K t {mt) + a e I e (m1?)I e (mt)\ {f < t < t ) (4.29) 

£=-oo 

where to = (e r °/4)ro. The first term in the sum corresponds to an expansion[6] of Ko(m\t — 1'|). 
The second term corresponds to the "reflected wave". The coefficients ag are to be determined 
by the continuity and boundary conditions. Similarly, 

oo 

G_+(r,r') = £ [h(rnt')K e+l (mt)-a e I f Xmt^ +1 (mt)}e^ +1 ^ (f < t < t ) (4.30) 

e=-oo 
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The corresponding elements of G are given by (4.27). There are similar expansions in the case 
t<t'< t . 

The coefficients ai will now be determined. When t' < t < t, (4.26) reduces to flG(r, r') = 
which means that G ++ is an analytic function of z and G_ + an antianalytic function. Therefore, 
as a function of z = re lip , G ++ is of the form 

oo 

G ++ = b e r e e ie< ? (t' < t Q < t) (4.31) 

e=-oo 

Taking into account the boundary condition G h = e llf G ++ at r = 2a gives 

G_+ = Y, b ^ 2a ) -^IW- ( t '< t °< t ) ( 4 - 32 ) 
e=-oo 

For a large disk, vq = 2atanh(ro/2) ~ 2aexp(— 2e _r °). The continuity of G ++ and G |_ at 

r = ro determines the coefficients ai and bf. One finds 

= exp[-(2l + l)2e- 7 P]^ + i(mt ) - /Q(mt ) 
° £ exp[-(2£ + l)2e- T -'']/ m (mt ) + /£(mt ) 1 ' j 

In the present large- fugacity limit, the Bessel functions in (4.29), (4.30), and (4.33) can be 
replaced by their asymptotic forms h(x) ~ (2nx)~ 1 / 2 e x and Kg(x) ~ (ir /2x) 1 / 2 e~ x . Then, 
whatever the relative order of t and t' might be, 

1 oo 

G ++ (r,r / )~- V { e - m l*-*'l -e- m ( 2t0 -*-* , Hanh[(2^+l)e- 7D ]}e i/v ' (t,t'<r ) (4.34) 

4vrt 

£=— oo 

and 

1 oo 

G_+(r, r') ~ V { e - m l'-''l + e -m&o-t-*) tan h[(2^ + l^ ]}^ 1 ^ (t, t' < r ) (4.35) 

47Tt 

£=— oo 

It should be recalled that these expressions are valid only near the boundary circle. 



4.3 Charge Fluctuations 

For the present model, by symmetry < Q >= and the variance of the total charge is 

(Q 2 )= f P ( V(r,r')dSdS' + f n(r)dS (4.36) 

Jr,r'<ro Jr<ro 

where p^(r, r') is the two-body charge density, n(r) the total particle density, and dS an 
area element on the pseudosphere. In the bulk, perfect screening is expected, and furthermore 
p( 2 \r, r') has a range in the geodesic distance between r and r' of the order of m~ l only. 
Therefore, the only contributions to (4.36) come from r and r' close to ro. When both r and r' 
are in the bulk (i.e. smaller enough than r ), p^ becomes a function p^ lk and it is convenient 
to define a surface part by p( 2 \r, r') = p^ k (r, r') + p^ rf (r, r'). Similarly, the density can be 
decomposed as n(r) = rabulk + w sur f(r). Assuming that perfect screening occurs in the bulk 
means 

J £ k (r,r')d5 + n bulk = (4.37) 
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where the integral extends on the whole pseudosphere. Using (4.37) allows to rewrite (4.36) as 
(Q 2 ) = ~ I &(r, r')dS dS' + [ £ f (r, r') + f n smi (r)dS (4.38) 

Because of the symmetry between positive and negative particles, p^(r, r') = 2[[/ + +(r,r / ) — 

U \-(r, r')] with the Ursell functions given by (4.25) and n(r) = 2n + {r) with n + given by (4.23). 

Thus, (4.36) becomes 

(Q 2 ) = -2m 2 I [|G ++ (r,r , )| 2 + |G_+(r,r')| 2 ]d5d5' + 2m I G ++ (r,r)dS (4.39) 

Jr,r'<ro J r<ro 

Using the Fourier series (4.34) and (4.35) in (4.39), and taking into account that in |G SS /| 2 only 
the term independent of ip survives the angular integration, gives 



Am 2 r 00 

(Q2 )= _ / y {e -2H*-*'l + e- 2m ^- t - t ')tanh 2 [(2£+l)e- T «]}d5d5 / 

+ — — / V {l-e- 2m (* -*Hanh[(2£ + l)e- T0 ]}d5 (4.40) 
4vrt Jt<tn „_ „ 



l=-oo 



(2) 

The first term in each sum corresponds to p^ k and nbulk, respectively, and therefore the second 
term corresponds to p^ r{ and n sur f, respectively. Using (4.38) rather than (4.36) gives instead 
of (4.40) 

(Q2 ) = _ y { _ f e - 2m ^dSdS' 

+ / e- 2m(2 * -*-*' ) tanh 2 [(2^+l)e- r "]d < Sd5'} 
Jt,t'<t 

oo „ 

y / e - 2m (*°-'Hanh[(2^ + I)e- T0 ]d5 (4.41) 

. . J t<tn 



2m 

£=-oo' 



47Tt ,— .Jt<to 



The integrands are indeed localized near the boundary circle. Using dS ~ to dt dip and perform- 
ing the integrations gives 

1 oo 1 oo 

= 4 E i 1 - tauh 2 [(2^ + l)e-^]} - - £ tanh[(2^ + l)e~^} (4.42) 

£=— oo £=— oo 

Finally, as tq becomes large, the sums can be expressed as integrals on the variable x = (21 + 
l)e _r °. Since tanhx is an odd function, the second sum can be considered as vanishing (actually, 
there are convergence factors at I — > ±oo, which have been omitted when the Bessel functions 
have been replaced by their asymptotic forms at fixed £). One is left with 

p-ro r°° p To 
(Q 2 ) = — (1 - tanh 2 x)dx = — (4.43) 
° J-oo 4 

in agreement with the macroscopic result (3.5), since here (5 = 2. 
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4.4 Surface Charge Correlations 



The first term in (4.34) or (4.35) corresponds to the bulk contribution (m/27r)Ko (m|t — t'|) or 
(m/2ir)e 1 ^ K\ (m \ t — t'|), respectively. The range mT x of these bulk contributions goes to zero 
in the large-fugacity limit. Thus, for t ^ t', only the second term survives. Let us assume that 
the relevant values of \ip\ are small. Since e~ r ° is small for a large disk, after (21 + l)e _T() has 
been replaced by 2£e~ T ° ~ eo£, the sum on t can be expressed in terms of an integral: 

tanh(e ^)e i£v ~ i / tanh(e ^) sm(£ip)d£ (4.44) 

Here too, there are convergence factors as £ — > ±oo, which have been omitted when the Bessel 
functions were replaced by their asymptotic forms at fixed £. These convergence factors can be 
taken into account by replacing tanh(eo^) by sinh(eo^)/ cosh(e^) (with e > eo), performing the 
integral which is a tabulated one[13], and taking the limit e — > eo afterwards. The result defines 
the integral as 

f 00 ivr 

i / tanh(eo^) sm(£cp)d£ = (4.45) 

7-oo eosmh^f 

The range in cp of this function is indeed of the order of eo, an a posteriori justification of the 
above assumption that \cp\ is small. Using (4.44) and (4.45), with | sinh(7r</?/2eo)| replaced by 
(1/2) exp7r|</?|/2eo), in (4.34) and (4.35) gives for the two-body charge density near the disk 
boundary 

p(2)( r , r ') = -2m 2 [|G ++ (r,r')| 2 + |G_+(r, r')| 2 ] = -^ e ^rn(2t -t-t') exp (44g) 

where eo ~ 2e _r ° and to ~ ae T °/2 have been used. This two-body charge density is indeed 
localized near the disk boundary. The surface charge correlation is defined as 

(a(if)a(O)} = f° dt f° dt'p (2) (r,r') (4.47) 

Using (4.46) in (4.47) and performing the integrals reproduces the macroscopic result (3.13), 
since here (3 = 2 and < a(ip) >= 0. 



5 ONE-COMPONENT PLASMA ON A PSEUDOSPHERE 

The macroscopic results (3.5) and (3.13) will now be checked on another solvable model, the 
one-component plasma. This is a system of one species of particles, of charges +1, embedded in a 
uniform background carrying the negative charge density — n&. At the inverse temperature (3 = 2, 
the system is exactly solvable in a variety of geometries, in particular for a large disk of radius 
R = aTQ on a pseudosphere [8]. A grand canonical ensemble is used. For the grand partition 
function to be convergent, it is necessary to define it with a fixed value [2] of the background 
charge density — n&; the fugacity £ controls the number of particles. Thus, in general, the system 
is not globally neutral, except for a particular choice of the fugacity. 

In the bulk the properties of the system are controlled by the background: the particle 
number density away from the boundary is rif,. However, near the boundary, the particle density 
differs from n&, and, since on a pseudosphere the neighborhood of the boundary has an area of the 
same order of magnitude as the whole area, this neighborhood gives an important contribution 
to the total number of particles and thus to the total charge of the system. 
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The macroscopic results (3.5) and (3.13) are expected to be valid only when the microscopic 
thickness of the surface charge density goes to zero. How to reach this regime in the most general 
way by varying both parameters and £ has not been clear to us. Here we content ourselves 
by considering the limit £ — > oo for a fixed value of n^. In this limit, the total charge of the 
system is expected to become infinite and to be carried by an infinitely thin surface layer. 

5.1 Summary of previous results [8] 

Again, for two points in the Poincare disk at r = (r, (p) and r' = (r', 0), one defines an auxiliary 
quantity G(r, r'), which now is just a scalar (instead of a 2 x 2 matrix). In the case of a large 
disk of radius R = aro, To — > oo, 

/ e ro+l \ a f r 2 \ ("+ 1 )/ 2 / r /2 \ (a+l)/2 ~ / rr , si e ify 

G(r ' r0 = C {—) V V~^J h W J l + WCe^ 

(5.1) 

where a = Aim^a 2 , x = 4£e T ° , and T(a,x) is the incomplete Gamma function 



poo 

T{a,x) = / t^e^dt (5.2) 

J X 

The particle number density was found to be 

ra(r) = G(r,r) (5.3) 

By a similar calculation, one finds for the two-body density Ursell function 

<7(r,r') = -|G(r,r')| 2 (5.4) 

In (5.1), r/(2a) = tanh(r/2). Only the case r large (r close to the boundary of the disk) 
will be needed. Then 1 - [r/(2a)f ~ 4e" T and [r/2a] e ~ exp(-2fe- T ). Let us assume that the 
relevant values oi\ip\ are small compared to 1. Then the sum on t can be replaced by an integral 
on x = 4^e _T ° . This gives for the density as a function of the distance (in units of a) from the 
boundary \ = tq — t 

n(A) = G(r,r)=Ce a e^W * ^— y (5.5) 



x u 



Integrating n(A) gives the average number of particles 

T(a, x) dx 



poo 

(N)=A(e a — 
Jo x 



+ 4na 2 (e a T(a, x) 



(5.6) 



5.2 Charge fluctuations 



For the one-component plasma the charge fluctuations are identical (for particles of charge 
+1) to the particle number fluctuations (Q 2 ) T = (N 2 ) T , since the background charge does not 
fluctuate. The charge fluctuations can be obtained either by integrating the correlation function 
(see eq. (5.4)) or by using the thermodynamic relation 



<«*> T -<^ (5.7) 
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This gives for a large disk 

(0 2 ) T T ^T(a,x) dx (5 8) 

4 io (x a + gT(a,x)) 2 

where we have defined the dimensionless parameter g = 47ra 2 Ce Q . For any finite value of ( 
and rib the integral in the last equation is different from 1, thus the predictions of macroscopic 
electrostatics are not satisfied. This is indeed expected since in general we are out of the validity 
domain of macroscopic electrostatics. As explained above, we expect the results from macro- 
scopic electrostatics to be valid only if the thickness of the layer of charge near the boundary 
is negligible compared to the macroscopic lengths: the radius of the disk R and the radius of 
curvature a. For the two-component plasma the thickness T of this layer is of order of the 
inverse of the fugacity vrT x = (AiraO^ 1 . For the one-component plasma we shall show that the 
situation is somehow different. 

Thus, before proceeding to study the charge fluctuations in the large-fugacity limit, let us 
study first how the thickness T of the charged layer near the boundary depends on g in this 
limit, since the situation is not as simple as it is for the two-component plasma. We will show 
that indeed T vanishes when g — > oo. 

For simplicity let us consider the case when a = 1. In units of a the thickness of the charged 
layer can be defined as the first moment of the density profile properly normalized 

^_ / "n(A)Ae~*dA _ J "n(A)Ae-*dA 

f °°n(A)e-MA n ^ 

where n = (N) j A is the average particle density. The e~ A factor comes from the area element 
dS = 27ra 2 sinhrdr near the boundary. For a = 1 the density profile (5.5) becomes 

f'°° nTP~ xeX 

n(X)=n b e 2X 9 ^ - dx (5.10) 
J x + ge x 

and the average density is given by 

nb Jo x + ge x 

Let us define x m as the principal solution of g = x rn e Xm ; incidentally, the function x m (g) is the 
Lambert function, which has many applications [15]. Now we write 

In the first integral (x < x m ) the second term in the denominator is negligible when g — > oo and 
then the integrand is 1. After the change of variable x — > x m + x, the second integral (x > x m ) 
is easily shown to have the limit In 2. This gives in the limit g ^ oo 

n 

— ~x m + ln2~x m (5.13) 
n b 

On the other hand, replacing expression (5.10) for the one-body density n(A) into the first 
moment of the density and performing the integral over A gives 



r n (X)Xe^dX = n b Tf^d* (5.14) 
Jo Jo x t 9 e 
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Again it is convenient to cut the integral in two intervals for x < x m and x > x rn . As in the 
case for n when g — > oo the second integral is negligible compared to the first. Then 



/ n(A)Ae~ A dA ~ n b / - 
Jo Jo 1 



e*r(o,_*)_ dx 



+ - 



n b 



fX m 

/ e x r(0,x)dx = 7 + e Xm r(0,x m )+lnx m (5.15) 
Jo 



where 7 ~ 0.577 is the Euler constant. Since, when x m — > 00, r(0,x m ) ~ e Xm /x m , the 
dominant term for the first moment of the density is the third term in the preceding equation 

roo 

/ n(A)Ae~ A dA ~ nf,lnx m (5.16) 
Jo 

Finally the thickness of the layer of charge near the boundary behaves as 

T~— lnx m ^0 (5.17) 

when x rn — > 00. Remembering that x m e Xm = g = 47ra 2 Ce Q one can notice that the dependence 
of the thickness on the fugacity is not trivial. It vanishes when ( — > 00 but very slowly contrarily 
to the case of the two-component plasma where T ~ 

Now we will proceed to prove that, in the limit g> 1, the charge variance (5.8) is equal to 
the prediction of macroscopic electrostatics (3.5). To be as general as possible we consider again 
any value of a. One can easily prove that the integrand in Eq. (5.8) is maximum when x = x m 
where x m is now given by g = x m /T(a, x m ) for any value of a. Doing the change of variable 
x — > x — x m in the integral (5.8) and replacing g by its expression in term of x m gives 

X-\-X m ) 



(QT = e i I ~) w (518) 

\ x m J T(a,x m ) 

When g^oowe have x m — > 00, T(a,x m ) ~ x^ _1 e _a;m and 



r( "' X + Xm) ^l + ^TV* (5.19) 



r(o;, x m ) y x 



Then 

(Q 2 ) T —± X _^]_A ^ _ dx (5.20) 

We notice that for large values of \x\ the integrand vanishes exponentially as e - ^'. Then, since 
x m — > 00, we can replace the lower limit of the integral by —00 and neglect x/x m in front of 1. 
This gives 

4 7^ (l + e-) 2 4 

Since = 2, this is the expected result (3.5) obtained from macroscopic electrostatics consider- 
ations. 
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5.3 Surface Charge Correlations 

Under the assumption that the relevant values of \p>\ arc small compared to 1, the same manip- 
ulations as the ones leading to (5.5) give 

G(r,r') = Ce Q e( Q+1 )^ / dx- (5.22) 

where g = 47ra 2 Ce a , A = tq — r, and A' = to — r'. 

Let us consider the case ip > 0. We are interested in the behavior of (5.22) as g — ► oo. This 
behavior will be shown to be determined by the pole of the integrand closest to the real axis 
in the upper half-plane. Let us assume that this pole has a large real part. Then, at this pole, 
T(a,x) behaves as the first term x a ~ 1 e~ x of its asymptotic expansion [13] and the denominator 
in the integrand becomes 1 + ge~ x jx. Let us look for a zero of this denominator at x = iir + x m . 
The equation for x m is (in + x m )e Xm = g, which becomes, in the large-g limit, x m e Xm = g. Thus, 
x m is large and real and the assumption that the pole has a large real part is a posteriori verified. 
The same reasoning gives other poles at x = (2n + l)i7r + x m , n £ Z. The residue of the pole at 
x = m + x m is easily found to have a modulus behaving as exp[— x m (e x + e A )/2] exp[— 7re T °ip/4]. 

Let I be the integral in (5.22). I is a part of an integral in the complex plane, along a contour 
C: C follows the positive real axis from to +oo, a large quarter of circle from +oo to +ioo, 
and comes back to the origin along the imaginary axis from +ioo to 0. The contribution from 
the large quarter of circle at infinity is easily seen to vanish, and therefore the contour integral 
is I — I' with I 1 the integral with x = iy pure imaginary varying from to ioo: 



roc - ly 2-^° yty 



r<=~ e •» 2 e y 4 
T = i / dy- (5.23) 



For large g, I' is easily seen to be of order 1/^7 — e Xm / x m . 

The theorem of residues says that I = I' + 2tt\ x sum of the residues of the poles inside C. 
I' is negligible (by a factor l/x m ) compared to the residue of the pole at iir + x m . The residues 
of the other poles have a factor exp[— (2n + l)7re T( V/4], n > 1 which makes them also negligible. 
A similar reasoning holds in the case ip < 0, and finally 

|G(r,r / )|~2vrCe Q e( a+1 )— exp[-x m ^— ]exp(-vr— M) (5.24) 

This form of |G| a posteriori justifies the assumption that the relevant values of ip are small 
compared to 1. Furthermore, in view of the fast decrease of |G| as a function of A or A' 
with a characteristic length l/x m (compare with the thickness of n(r) which was found to be 
(l/x m )lnx m ), a simpler form is 

A A' e T ° 
|G(r,r')| ~ 27rCe Q exp[-x m (l + - + — )] exp(— tt— |^|) (5.25) 

The Ursell function is obtained by using (5.25) in (5.4): 

e T ° 

U(r, r') ~ -(2<e a ) 2 exp[-x m (A + A')] exp(-2x m ) exp(-7r— \ip\) (5.26) 
The surface charge correlation is defined as 

(ct(v?)ct(0)) t = a 2 dX dA'[/(r,r') (5.27) 
Jo Jo 
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Performing the integrations and using 4ira 2 (e a = x m e Xm reproduces the macroscopic result 
(3.13) at j3 = 2. 

6 CONCLUSION 

The charge fluctuations for a two-dimensional classical Coulomb fluid are drastically changed 
by the introduction of a negative curvature of space. 

In the case of a flat disk communicating with a reservoir (grand-canonical ensemble), the 
total charge Q essentially does not fluctuate (bringing an additional charged particle from infinity 
would cost an infinite energy). In the macroscopic limit, one can define a surface charge density 
a (charge per unit length on the boundary circle). The two-point correlation function of a has 
an algebraic only decay (1.1), behaving as the inverse square distance between the two points 
(while the charge correlation function in the bulk has a faster than algebraic decay). 

In the case of a disk on a pseudosphere (an infinite surface of constant negative curvature), 
in the macroscopic limit, the total charge Q does fluctuate with the variance (3.4). Furthermore 
the two-point correlation function of the surface charge density a has a fast (exponential) decay 
(3.13) as a function of the angular distance \ip\ between the two points. 

This change of behavior of the surface charge correlation is related to the well-known fact 
that a negative curvature acts as a mass in the field equations. The curvature replaces the flat 
logarithmic Coulomb potential by the potential (2.11) which has an exponential decay at large 
distance s. For a flat disk, the algebraic decay of the two-point surface charge correlation is 
due to these field lines which connect the two points through the vacuum outside the disk. On 
a pseudosphere, these field lines outside the disk nevertheless carry an exponentially decaying 
interaction. 

For retrieving the macroscopic limit from microscopic models, it is necessary that the thick- 
ness T of the surface charge density be negligible compared to the macroscopic length scales. On 
a pseudosphere with a radius of curvature a, in a disk of radius R, we have considered only the 
case iJ> a. The macroscopic behavior is expected to hold only when o>T. The two exactly 
solvable microscopic models which have been considered do exhibit the expected macroscopic 
features when this condition is satisfied. 

ACKNOWLEDGEMENT 

J.M.Caillol brought to our attention the literature about the Lambert function. The au- 
thors acknowledge support from ECOS-Nord/COLCIENCIAS-ICFES-ICETEX action C00P02 
of French and Colombian cooperation. 

References 

[1] E.H.Lieb and J.L.Lebowitz, Adv. Math. 9:316 (1972). 

[2] E.H.Lieb and H.Narnhofer, J.Stat.Phys. 12:291 (1975). 

[3] B.Jancovici, J.Stat.Phys. 110:879 (2003). 

[4] B.Jancovici, J.Stat.Phys. 80:445 (1995). 

[5] B.Jancovici, J. Phys.:Condens. Matter 14:9121 (2002). 



18 



[6] A.Erdelyi et al, Higher Transcendental Functions (McGraw-Hill, New York, 1953), Vol.11. 

[7] A.Erdelyi et al, Tables of Integral Transforms (McGraw-Hill, New York, 1954), Vol.I,p.387. 

[8] R.Fantoni, B.Jancovici, and G.Tellez, J.Stat.Phys. 112:27 (2003). 

[9] B.Jancovici and G.Tellez, J.Stat.Phys. 91:953 (1998). 
[10] P.J.Forrester, J.Chem.Phys. 95:4545 (1991). 
[11] B.Jancovici and G.Tellez, J.Stat.Phys. 82:609 (1996). 
[12] F.Cornu and B.Jancovici, J.Chem.Phys. 90:2444 (1989). 

[13] S.Gradshteyn and I.M.Ryzhik, Table of Integrals, Series, and Products (Academic, New 
York, 1965). 

[14] J. Zinn- Justin, Quantum Field Theory and Critical Phenomena (Clarendon, Oxford, 1989). 

[15] R.M.Corless, G.H.Gonnet, D. E.G. Hare, D.J.Jeffrey, and D.E.Knuth, Adv.Comput.Math. 
5:329 (1996). 



19 



